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ABSTRACT 

We investigate a number of potential foregrounds for an ambitious goal of future radio 
telescopes such as the Square Kilometer Array (SKA) and Low Frequency Array (LO- 
FAR): spatial tomography of neutral gas at high redshift in 21cm emission. While the 
expected temperature fluctuations due to unresolved radio point sources is highly un- 
certain, we point out that free-free emission from the ionizing halos that reionized the 
universe should define a minimal bound. This emission is likely to swamp the expected 
brightness temperature fiuctuations, making proposed detections of the angular patch- 
work of 21cm emission across the sky unlikely to be viable. An alternative approach 
is to discern the topology of reionization from spectral features due to 21cm emission 
along a pencil-beam slice. This requires tight control of the frequency-dependence of 
the beam in order to prevent foreground sources from contributing excessive variance. 
We also investigate potential contamination by galactic and extragalactic radio re- 
combination lines (RRLs). These are unlikely to be show-stoppers, although little is 
known about the distribution of RRLs away from the Galactic plane. The mini-halo 
emission signal is always less than that of the IGM, making mini-halos unlikely to be 
detectable. If they are seen, it will be only in the very earliest stages of structure for- 
mation at high redshift, when the spin temperature of the IGM has not yet decoupled 
from the CMB. 

Key words: cosmology:theory - galaxies:formation - large-scale structure of universe 



1 INTRODUCTION 

We have as yet no direct observational probes of neutral gas 
beyond the epoch of reionization. One promising technique 
which has attracted much attention is 21cm tomography of 
neutral hydrogen (Scott & Rees 1990; Madau, Meiskin & 
Rees 1997; Tozzi et al 2000). 21cm emission and/or absorp- 
tion from neutral hydrogen at high redshift should exhibit 
angular fluctuations as well as structure in redshift space. 
These fluctuations are due to spatial variations in the hy- 
drogen density, ionization fraction, and spin temperature, 
and may be detectable by future radio telescopes such as the 
Square Kilometer Array (SKA)^ and the Low Frequency Ar- 
ray (LOFAR)^. Although the energy density in 21cm emis- 
sion is about two orders of magnitude less than the cosmic 
microwave background (CMB), the extreme smoothness of 
the CMB both spatially and in frequency space would allow 



^ see 'http:/ /www. nfra.nl/skail 
^ see http://www.astron .nl/lofarl 



this signal to be teased out. The expected brightness tem- 
perature fluctuations on arminute scales is about two orders 
of magnitude larger than CMB fluctuations. In principle, 
21cm tomography would allow us to map out the topology 
of reionization. This could indirectly constrain the nature of 
the ionizing sources, telling us whether these sources were 
faint and numerous or bright and rare. It would thus be 
an invaluable tool for probing the state of the intergalactic 
medium during an epoch when traditional probes, such as 
Lya absorption, fail (Gunn-Peterson absorption fully satu- 
rates at hydrogen neutral fractions xhi ~ 10"''). 

We present a study of possible foregrounds for these ob- 
servations. On the arcminute angular scales where the signal 
is expected to peak, detailed studies of the CMB (involving 
extrapolation from higher frequencies) tell us that fluctu- 
ations in the galactic foreground emission are likely to be 
relatively unimportant. These galactic foregrounds are fairly 
well understood and multi-frequency observations should en- 
able us to subtract out both the free-free and synchrotron 
components (see Shaver et al (1999) for a detailed discus- 
sion). Unresolved extra-galactic radio sources, on the other 
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hand, could give rise to brightness temperature fluctuations 
which would swamp the expected signal, as discussed by Di 
Matteo et al (2002). However, this model for source counts 
was based on surveys with limiting flux densities ~ lOOmJy 
at 150 MHz, which were then extrapolated down five or- 
ders of magnitude past ~ l/iJy (the expected limiting point 
source sensitivity of SKA). Therefore, as they acknowledge, 
this estimate of foreground contamination is highly uncer- 
tain; a turnover in source counts at lower flux levels could 
significantly reduce foreground contamination. 

We present a minimal estimate of the brightness tem- 
perature fluctuations, based on free-free emission from the 
ionizing sources that reionized the universe (Oh 1999). We 
show that brightness temperature fluctuations from these 
sources exceed the expected signal, making it unlikely that 
the previously proposed angular 21cm tomography is fea- 
sible. We shall show that this estimate is relatively model- 
independent and depends primarily on the integrated ioniz- 
ing emissivity. 

An inability to detect angular brightness temperature 
fluctuations need not render 21cm tomography studies pow- 
erless. Since the foreground signal is expected to be smooth 
in frequency space, 21cm spectral features along a pencil- 
beam slice, corresponding to alternating patches of neutral 
and ionized hydrogen, could yield invaluable information on 
the topology of reionization. However, this relies heavily on 
the lack of foreground spectral contaminants. We investi- 
gate two possible contaminants: spectral fluctuations in the 
foreground due to the frequency- dependent beam, and radio 
recombination lines from galactic and extragalactic ionized 
gas. The former can probably be dealt with if the frequency 
dependence of the beam can be controlled and its side-lobes 
mapped out. The latter is much more uncertain but is un- 
likely to be a show-stopper. If, in fact, it does turn out to 
be important, it will yield the unexpected bonus of new in- 
formation about the clumping of ionized gas in our galaxy 
and in the universe as a whole. 

In all numerical estimates, we assume a 
ACDM cosmology where (f2m, Ha, fibft^, /i, ag^^-i) = 
(0.3,0.7,0.019,0.7,0.9). 



2 FREE-FREE EMISSION FROM IONIZING 
SOURCES AS A FOREGROUND 

We begin by summarizing the properties of the expected 
signal; the reader is referred to the literature (Scott & Rees 
1990; Madau, Meiskin & Rees 1997; Tozzi et al 2000) for de- 
tails. Consider a patch of IGM, with spin temperature Ts, 
which fully fills the beam and has a radial velocity width 
larger than the bandwidth of the radio telescope. The dif- 
ferential brightness temperature between this patch and the 
CMB is (Tozzi et al 2000): 

-~'-(S)(^)'"(^)^« 

Of course, if ionized bubbles exist within this patch then the 
signal will be diluted by the corresponding filling factor. 



Initially, Ts = Tcmb, but over time Lya photons from 
the soft UV background emitted by the first stars and 
quasars will couple the spin temperature of the gas to its 
kinetic temperature through the Wouthuysen-Field effect. 
Since initially Tk < Tcmb, neutral hydrogen will be visible 
for a brief period in absorption until the same Lya photons 
also heat the gas through recoil of the scattered Lya photons 
and make Ts > Tcmb - At this point the signal is only visible 
in emission; from equation Q the brightness temperature is 
roughly independent of the spin temperature. Hereafter, we 
shall focus solely on the emission signal. 

If linear theory is used to compute expected density 
fluctuations, the brightness temperature fluctuations peak 
on arcminute scales and have a rms value of {ST^y/^ 
lOmK for a ACDM cosmology (see Fig 1 of Tozzi et al 
(2000)). To give a sense of the lengthscales involved, a 
bandwidth Av corresponds to a comoving length L ~ 

(1 + z)cH{zyAu/i' = 8.6 m^y'^^^ (A///lMHz)/i-i Mpc, 
where v = 1-4/(1 -|- 2)Ghz is the observation frequency, and 
an angular diameter AQ corresponds to a comoving trans- 
verse length L = A6i(l + 2)/dA(z) ~ 1.9 (AS/l') /i"^ Mpc at 
z=9. 

The contribution of radio-loud AGN and radio galax- 
ies at low flux levels is highly uncertain. We therefore con- 
struct a minimal model of the low-frequency radio back- 
ground from free-free emission by ionizing sources. This 
background is unavoidable in the sense that a minimal emis- 
sivity in ionizing photons is necessary to reionize the uni- 
verse. We present a brief summary here; refer to Oh (1999) 
for details. The free-free luminosity of an ionizing source 
can be estimated from the fact that the free-free emissiv- 
ity e^/ oc rig is directly proportional to the recombination 
rate and thus to the production rate of ionizing photons: 
Ll^ oc {nl)V oc Ar'-='=°™'' oc (1 - f,^sc)Nion. We thus obtain: 

» "^'»" ( io..p£;:„.- ) (H^) 

where Nion is the production rate of ionizing photons, and 
fesc is the escape fraction of ionizing photons from the host 
galaxy. For a Salpeter IMF with stars of solar metallicity, 
Nion = lO'^^photonss"^ (SFR/lMoyr"^) . For a given star 
formation rate (SFR), stars with zero metallicity can be 
one to two orders of magnitude more effective at producing 
ionizing photons, since the effective temperature of these 
stars is higher (Tumlinson & ShuU 2000), and the IMF 
of zero-metallicity stars is often thought to be top-heavy. 
The escape fraction of ionizing photons is observed to be 
fesc ~ 5% in the local universe (Leitherer et al 1995), and 
thought to decline with redshift (Wood & Loeb 2000; Ri- 
cotti & ShuU 2000), although this is highly uncertain (for 
calculations arriving at the opposite conclusion, see Fujita 
et al (2002)). The observed free-free ffux is then given by 
Sa = L^{l + z)/4ndl,or. 

" ( ios3p^;:uss-0 (^) ""''-^'^ 

Note that the uncertainty in the escape fraction of ionizing 
photons introduces uncertainties of at most a factor of only 
a few in the free-free luminosity. Only in the unlikely sce- 
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nario of virtually all the ionizing photons escaping into the 
IGM without any photo-electric absorption in the host ISM, 
would order of magnitude uncertainties creep in as feac — > 1. 

To compute the brightness temperature fluctuations 
due to these sources, we need to construct a luminosity 
function. Oh (1999) followed the star formation model 
of Haiman & Loeb (1997), who assumed that some con- 
stant fraction of the gas /, in halos with virial tempera- 
tures Tvir > lO'^K fragments to form stars, and that each 
starburst lasts for ~ lO'^yr. The star formation efficiency 
parameter /* is tuned for consistency with the observed 
metallicity of the IGM at z=3, with /* ~ 1.7 - 17% cor- 
responding to Z ~ 10"'' — W~^Zq, respectively. With this 
prescription, Nion « 2 x 10^^ [Mhaio/W"^ Mq) photonss"\ 
Press-Schechter theory can be used to compute the abun- 
dance of halos and hence the source counts dN/dS. The 
power spectrum of Poisson fluctuations can then be com- 
puted from: 



C 



' ,c:d.N 2 



(4) 



where Sc is the minimal flux above which point sources can 
be identified and removed. The power spectrum due to clus- 
tering can be computed from: 



c; 



clustering 
I 



r2 



(5) 



where wi is the Legendre transform of the angular correla- 
tion function of sources w{9), and 7^ = J^" dS{dN/dS)S is 
the surface brightness of the radio background at frequency 
f. The spatial correlation length can be computed from lin- 
ear theory assuming linear bias; the decrease in the linear 
growth factor at high-redshift is increased by the number- 
weighted bias of objects with T^ir > 10*. This results in 
an angular correlation function w{6) = {0/0o)~^'^ with a 
nearly constant angular correlation length 9o ~ 2' indepen- 
dent of redshift for z > 3 (see Fig. 5 of Oh (1999)). The 
rms temperature fluctuations can then be computed from 
Trms = + l)C;/47r)^''^ X (? l2kBv'^, where the frequency 
dependence of the temperature conversion is appropriate in 
the Rayleigh- Jeans limit. This quadratic frequency depen- 
dence implies that the population of free-free emitters only 
becomes important at low frequencies. Thus, for instance, 
their existence does not violate any present CMB distortion 
constraints. 

What is the minimal flux above which point sources 
can be identifled and removed? The detector noise is: (Rohlfs 
& Wilson 1996): 



'5'inst — 



2fcre, 



/ Ai/ \ "^/^ / t \ -1/2 



AftV2tA 

where we have used Aeft/Tays = 2 x lO^cm^K^^. On the 
other hand, the confusion noise is: 



5'conf=(Cr'" = **=''(Se)e^)'^'^2^ ^' 



20nJy 



03^) 



(1999)). This is because Cf'""*'='' cx /^(Sc), and from equa- 
tion the majority of sources which constitute the free-free 
background are generally fainter than ~ 20 njy; removal of 
rare bright sources has little effect on the mean free-free 
background ly. We have assumed unresolved point sources 
so that Q is given by the beam resolution; the confusion noise 
becomes correspondingly larger if the source is extended. 

Confusion noise, rather than instrumental noise, is 
the limiting factor in source identiflcation and removal for 
free- free sources. The free- free flux is largely frequency- 
independent. If instrumental noise were the limiting factor, 
point sources could be identifled and removed with higher 
frequency observations, for which instrumental noise is much 
lower (for instance, in 10 days a 5cr detection of a 16 nJy 
source at 2 GHz is possible, for a bandwidth IS.V = 1 GHz). 
However, since confusion noise is frequency independent, 
once the confusion limit is reached no further source removal 
is possible, even with multi-frequency observations. 

We compute the Poisson and clustering contributions 
to the free-free background with the above model. The re- 
sults are displayed in Figure Q The 21cm signal is likely to 
be swamped by Poisson fluctuations at small scales and by 
fluctuations due to clustering of free- free sources at all scales. 
This calculation should be compared with that of Di Matteo 
et al (2002), who use a power law extrapolation of observed 
number counts to compute the same quantities. Although 
we arrive at similar conclusions, below we argue that these 
estimates (in particular the temperature fluctuations due to 
the clustering component) are significantly more robust. 

The results obviously depend on the specifics of the star 
formation model assumed. Since this is highly uncertain, we 
highlight here the most robust features of the calculation, 
which are largely model independent. The number counts in 
the Poisson signal are dominated by rare bright objects just 
below the detection threshold Sc., and can vary widely de- 
pending on the shape of the luminosity function at these low 
fiux levels. We therefore disregard this term as it is exces- 
sively model dependent. On the other hand, the clustering 
term, equation © is significantly more robust. Our estimate 
of Qo ~ 2' is likely a minimal estimate of the angular clus- 
tering strength since it is dominated by the bias of halos 
at the threshold virial temperature T^ir ~ lO^K, which are 
most abundant. Star formation in lower virial temperature 
halos, which cannot cool by atomic line cooling, is expected 
to be extremely inefficient; see e.g. Barkana & Loeb (2001) 
for a general review. If, in fact, star formation is more effi- 
cient in deeper potential wells (e.g., since feedback processes 
are less devasting there), the luminosity weighted correlation 
function will be more heavily weighted towards rare massive 
halos, which are more highly biased. Thus, the wi term is 
generally a minimal estimate. The surface brightness term 
it is independent of the model for the luminosity function 
and depends only on the overall normalization of the comov- 
ing luminosity density. We can see this from the equation of 
cosmological radiative transfer (Peebles 1993): 



Thus, for a lOcr detection, the cutoff flux is Sc ~ 20nJy. 

Note that Sconf oc (Cf''"'*'=''(S'c))''"^ cx S^'^ depends only 
weakly on the fiux cutoff (see Fig. Q and also Fig. 4 of Oh 



dz 



dt e[uo{l + z)] 
'dz {l + z)3 



(8) 
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Figure 1. Brightness temperature fluctuations induced by free- 
free emission from ionizing sources, compared with the expected 
21cm signal, at a frequency of 150 MHz (corresponding to an 
observed redshift of 2 = 8.3), and assuming a star formation effi- 
ciency corresponding to a Lyo forest metaUicity of Z = 10~^'^Zq 
at z = 3. SoUd and dashed lines assume point source removal 
down to 20 njy and 2/iJy respectively. The 21cm signal is likely 
to be swamped by fluctuations due to clustering of the free- 
free sources at all scales. The estimate for the clustering com- 
ponent is fairly robust and model independent (see text). Note 
that T^l 



(X n,(2 = 3) oc Z{z = 3). 



where e{i>) is the comoving emissivity. Since e oc Q,{1 — fesc) 
(from equation the radio surface brightness I{i^o) oc 
f2«(l — /esc)/(l + z)^, where z is the median redshift at 
which ~ 50% of all stars have formed. Thus, a given co- 
moving stellar density fi* directly and robustly implies a 
minimal free-free surface brightness on the sky (the depen- 
dence on the star formation history f2,(2) is weak, since 
most stars formed at late times). In our calculation we have 
conservatively included only sources at 2 > 3 (which are 
too faint to be identified and removed) and simply normal- 
ized to n^z = 3) = 6 X 10"^ (Z/IQ-^-^Zq) Qb, where Z 
is the mean metaUicity of the IGM at 2 ~ 3, and we as- 
sume a Salpeter IMF where ~ IMq of metals form for every 
~ IOOMq of stars formed. As previously noted, because the 
majority of sources are significantly fainter than the thresh- 
old flux for point source removal Sc, the mean surface bright- 
ness depends only weakly on Sc (see Fig. 4 of Oh (1999)). 
Since both the estimates for the clustering strength wi and 
the mean surface brightness I^, are minimal and robust, we 



conclude that our estimate for C? 



constitutes a robust 



minimal estimate for angular temperature fluctuations on 
the sky. These fluctuations will swamp the expected bright- 
ness temperature fluctuations due to 21cm emission by at 
least an order of magnitude. 



3 FOREGROUNDS FOR SPECTRAL 
MEASUREMENTS 

Even if brightness temperature fluctuations due to 21cm 
emission are not detectable, spectral variation along a 
pencil-beam might still be detected. This relies on the 
fact that foreground sources should in general have smooth 
power-law continuum spectra, and averaging over many 
sources with different spectral indices and spectral struc- 
ture will yield a smooth power-low foreground (see Shaver 
et al (1999) for simulations of this) which would still allow 
detection of spectral structure corresponding to alternating 
regions of neutral and ionized hydrogen. To give an idea of 
scale, an ionized 'bubble' created by a source forming stars 
at a rate ~ 2M0yr~^ for ~ lO^yr has a diameter ~ 1.4[(1 -I- 
z)/10]~^/i~^Mpc comoving, while a bandwidth Av spans 

a comoving length L ^ 8.6 (^r)"^''^ {Au/lMRz)h-'^ Mpc. 
Of course, as reionization proceeds the ionized regions ex- 
pand in size and neutral regions contract. As Di Matteo et 
al (2002) point out, observing with 2 MHz resolution at 
I'o = 150MHz with an error of A/3 ~ 0.05 in the measured 
foreground spectrum spectral slope should allow one to dis- 
tinguish a signal that is [{u/vo)^'^ — 1] ~ 6 x 10~* times 
smaller than the foreground. Here, we discuss two potential 
foregrounds for such spectral measurements: beam smearing 
of foreground sources and radio recombination lines. 

The radio telescope does not sample exactly the same 
patch of sky at all frequencies, since the beam-size and side- 
lobes are frequency dependent. This would by itself natu- 
rally introduce frequency space fluctuations in the measured 
brightness temperature as additional sources fill the beam, 
an effect of order: 



AO, AO Av 
~ 2 — ~ 2 

Q. e V 

1.4xl0-^f^ 
V2MHz 



(9) 



10 



Since we have 



seen that the clustering foreground 
{T'fgY^'^^lQ{T2xcm)^^'^j the noise in frequency space in- 



troduced by beam-size variation will be STj 



fa 



few X 



0.1{r2icm)^''^) comparable to temperature fluctuations due 
to 21cm emission. An exception may be in the very 
early/late stages of reionization, when the size of ion- 
ized/neutral patches is small; the beamsize does not vary 
signficantly over the small bandwidth Av necessary to de- 
tect spectral feaures. Otherwise, it would be necessary to 
develop scaled beams and controlled sidelobes, which are 
capable of sampling exactly the same patch of sky at all 
frequencies. Note that frequency calibration of the telescope 
for pencil beam tomography is considerably more challeng- 
ing than frequency calibration for detecting an all-sky signal, 
such as that produced at the tail end of reionization (Gnedin 
& Ostriker 1997). In the case of an all-sky signal, one could 
simply average the results over many independent patches 
of sky, as well as perform a differencing experiment off the 
moon, which blocks emission behind it (Shaver et al 1999). 

Another possible source of structure in frequency space 
is the signal from radio recombination lines (RRLs). These 
lines have been observed both in emission and absorption 
at a wide range of frequencies, although generally only 
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(10) 



against bright sources in the Galactic plane (see Gordon 
& Sorochenko (2002) for a comprehensive review). The fre- 
quencies of the hydrogen lines are: 

v^l53An(^) ^MHz. 
V350y 

The important question is whether RRLs will provide a sig- 
nificant source of contamination on lines of sight outside 
the Galactic plane. There have been few RRL line searches 
away from bright sources outside the Galactic plane, and 
in any case surveys for radio recombination lines generally 
have low-frcqucncy detection limits above the strength of 
the 21cm features we seek. We can place a useful limit using 
the fact that the observed brightness of both RRL and Ha 
emission depends on the emission measure EM = / dsn'i, 
and that optical surveys for Hq have much lower limiting 
sensitivities. Fabry-Perot surveys have detected Ha emission 
from every Galactic latitude, with minimal values 0.25 — 0.8 
Rayleighs toward the Galactic pole (Reynolds 1990). The 
emission measure corresponding to an observed Ha intensity 
Ic, (in Rayleighs) is; EM{Ha) = 2.75Ti'-3/c,cm"^pc. Based 
on the width of the Ha line, Reynolds (1990) estimates 
T ~ 8000K. Given the density, temperature and size of a re- 
gion, as well as an estimate of the linewidth and a frequency 
to observe at, one can also calculate the recombination line 
optical depth (Shaver 1975): 

Te dlnbr, 



for reasonable assumptions about the excitation parame- 
ter only nearby sources {D < lOMpc) will be detectable in 
RRL emission if spontaneous emission is responsible (Shaver 
1978). It was thought that stimulated emission could allow 
observation of much more distant sources (Shaver 1978), 
but that has not turned out to be the case. 

Extragalactic RRLs could be a foreground in the search 
for 21cm absorption against radio-loud sources at high red- 
shift (Carilli, Gnedin & Owen 2002; Furlanetto & Loeb 
2002). Since high-redshift objects are in general much denser 
rie oc (l+z)^, their emission measures are much higher. Con- 
sider an isothermal disk at temperature Tgas embedded in 
a halo with virial temperature Tvir- One can solve for the 
density profile of such a disk, n oc exp(— r/i?d)scch^(2/zo) 
(Wood & Loeb (2000); see also Oh & Haiman (2002)), 
and the emission measure of the disk will be EM{r) « 
f dzn^{r,z). For a disk seen face-on, the emission measure 
becomes a function of the radial distance r from the center: 



i.M(0 ^2x 10^ ( A) ^^(^)-^( 



T 

vir 



lo:o5j l^j ^"K-^j 



104 

2r 



5 X 104 K 



3/2 



(12) 



518EM , 

TL = — TT^On 



20.8- 



dn 



where fd = {Mdisk/Mhaio)/{^b/^o) is the fraction of 
baryons in the disk, A is the spin parameter and Rd ~ 
I\f2. The optical depth at line center for a An = 1 
RRL will be: 



where EM is the omission measure in cm~®pc defined by 
ri^ I, v is the line center frequency, Vl is the line width in 
km/s, Te is the electron temperature, and 6„ is the depar- 
ture coefficient, which takes into account departures from 
local thermodynamic equilibrium (LTE). For LTE, &„ = 1. 
For an unresolved line, the velocity width is simply: AVl = 
c{Ai'/u), where Av is the observational bandwidth. Assum- 
ing /„ = 0.25Ji (which corresponds to EM = 0.5cm~^pc) 
and Tl = tlTe, we can employ the optical depth calculation 
to estimate the hydrogen RRL temperature: 



TL = 2x10" 



-5/2 



EM 



108 



cm~^pc 



0( 



i^)"^<-) 



Tl = 1.3x10" 



EM 



assuming local thermodynamic equilibrium and that the line 
is unresolved. This could be considerably enhanced if stim- 
ulated emission and/or non-LTE effects become important 
(which is likely since one is looking along the line of sight to a 
radio- bright source). This optical depth is comparable to the 
central optical depths due to 21cm absorption by the IGM 
(Carilli, Gnedin & Owen 2002) or mini-halos/mini-disks 

If 



( sOOOg ) ' ( o5cti^6 c ) (l^lz) K(ll) (Furlanetto & Loeb 2002), which are of order r ~ 10 



which is negligible. Note, however, that this estimate as- 
sumes thermodynamic equilibrium. Particularly for lower 
frequency lines, stimulated emission may become important. 

Carbon radio recombination lines axe another possible 
source of contamination. They have been detected in the 
frequency range 34-325 MHz toward Cass A in the Galactic 
plane (Payne et al 1989), and also in a 327 MHz survey 
centered at 6 = 14° (Roshi et al 2002). Optical depths are 
generally r ~ 10~^ and they are thought to arise in dense, 
cold (Te ~ 20 — 200K), partially ionized regions where non- 
LTE effects are important. Carbon RRL emission is gener- 
ally confined to 6 < 3° and their intensity out of the plane 
is probably small (Roshi et al 2002). 

Detections of radio recombination lines from extragalac- 
tic sources have generally been confined to bright starburst 
galaxies and the intensity appears to be correlated with the 
star formation rate (Phookun et al 1998). They are un- 
likely to be significant contaminants. One can show that 



one uses Press-Schechter theory to calculate abundances of 
halos with T^ir > lO'^K, a random line of sight would inter- 
sect ~ few x 0.1 disks per redshift interval, for the redshifts of 
interest (Furlanetto & Loeb 2002) . Even a single disk along 
the line of sight would introduce a plethora of RRL features 
in emission and/or absorption. This would make the task 
of picking out the spectral features due to 21cm absorption 
much more difficult. On the other hand, such observations 
would also be an invaluable probe of gas clumping at high 
redshift. However, this is a rather speculative scenario, and 
we shall not comment on it further. 

In general, radio recombination lines are unlikely to be 
significant contaminants. If they are, they will have to be 
removed with observations of higher spectral resolution (us- 
ing the fact that they occur at known frequency intervals 
for identication and removal). In many ways their detection 
could be an unexpected bonus, as RRLs have the potential 
to teach us a great deal about the clumping of ionized gas 
both within our Galaxy and in the universe as a whole. 
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4 DETECTION OF MINI-HALOS 

Recently, it has been proposed that mini-halos with T < 
lO^K, which are not colhsionally ionized but which have 
gas densities sufficiently high to coUisionally couple the spin 
temperature to the kinetic temperature, may be detectable 
in emission (Iliev et al 2002a; Iliev et al 2002b). For Ts > 
TcMB, the flux is independent of the spin temperature and 
depends only on the HI mass. Although the signal from an 
individual mini-halo is very faint, the combined signal from 
many mini-halos within a sufficiently large comoving volume 
may be detectable. We point out that because Si, oc Mhi, 
the flux from mini-halos will always be swamped by the flux 
from the neutral IGM, as the HI mass in the IGM is larger. 
The ratio of fluxes is given by the relative fraction of HI in 
collapsed halos with Tvir < lO^K: 

^ ^ J_ dM^M (14) 

where dN/dM is the halo abundance given by Press- 
Schechter theory, Mmgh = 10**[(1 -I- 2)/lO]"^/^M0 is the 
mass of halos with T^ir = IC'K, and Miow is the baryonic 
Jeans mass below which gas cannot accrete onto halos. If 
there is no heating of the IGM, Miow is given by the cosmo- 
logical Jeans mass Mj = 8 X 10^ [(1 + z)/lO]^/^M0, where 
the temperature of the IGM Tigm = 1.7[(1 -I- 2)/10]^K is 
simply set by the amount of adiabatic cooling since decou- 
pling from the CMB. In reality it is likely that the gas will 
be heated up to higher temperatures by soft X-ray and Lya 
photon recoil heating from the first sources of light, as well 
as by free-free emission from collapsing structures (Madau, 
Meiskin & Rees 1997). This will suppress the accretion of 
gas onto the smallest mini-halos. 

We plot the ratio of fluxes from equation 11411 in Fig Q . 
Even if heating is negligible, this ratio is at most a few x 0.1. 
This may be amplified in high density peaks by a factor of 
no more than a few due to clustering bias (Iliev et al 2002b) . 
Thus, the flux from mini-halos only dominate in two limits: 
(i) in reionized regions, for self-shielded mini-halos which 
have not yet been photo-evaporated, and (ii) very early in 
the structure formation process, when the background ra- 
diation field is low and the spin temperature of the IGM 
is still coupled to that of the CMB. The former case is un- 
likely to be detectable as the mini-halos contribute a small 
signal superimposed on the much larger fluctuating signal 
of alternating neutral and ionized regions of the IGM. The 
latter may be detectable provided the mini-halos indeed re- 
main neutral and do not form stars via H2 cooling: even a 
single star in a halo is capable of photo-ionizing and photo- 
evaporating much of the gas. 

It is worth checking if this is possible. From Haiman, 
Abel & Rees (2000) (see their Fig. 6), the level of back- 
ground UV flux capable of photo-dissociating H2 in mini- 
halos and preventing star formation is: 

Jmss ~ 10"^ - 10"^ J21 (15) 

in the redshift range z=10-20 where J21 (in units of 
10~^^ergs~^cm~^Hz~^sr~^) is the UV flux in the Lyman- 
Werner bands 11.2— 13.6eV. On the other hand, the critical 
thermalization flux of Lya photons (which consists of red- 
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Figure 2. Relative 21cm flux from mini-halos compared to IGM 
in a neutral patch which fills the beam, assuming Tg ^ Tcmb 
for both. Since S^)'^'"" / Sl'^'^' = flhalos/^b, the 21cm emission 
from the IGM will always be greater except in the earliest stages 
of structure formation when the spin temperature of the IGM is 
still coupled to that of the CMB. The relative mass fraction (and 
hence 21cm flux) in mini-halos decreases still further if the IGM 
is heated to 'warm' temperatures (e.g. Tjqm ~ 200K, 2000K as 
shown here) by soft X-rays or atomic recoil from Lya photon scat- 
tering, which increases the baryonic Jeans mass and suppresses 
accretion onto mini-halos. 

shifted photons from the Lyman- Werner bands) required to 
decouple the spin temperature of the IGM from the CMB 
and drive Ts Tk is (Madau, Meiskin & Rees 1997): 

Jtherm ~ 0.9J21 [^^) ■ (16) 

Since Jdisa ^ Jtherm, there certainly exists a period in the 
history of the IGM where gas in mini-halos cannot cool to 
form stars, but the IGM does not emit in 21cm radiation. In 
this case, the halos would give rise to a period in which 21cm 
is detectable in emission before it is seen in absorption when 
Ts.iGM < Tcmb- 21 cm radiation would once again be de- 
tectable in emission due to heating when Ts.igm ^ Tcmb- 
The duration of each of these epochs can be used to time 
the rise of the radiation field. Of course, this is all assuming 
that the noise from foregrounds previously discussed can be 
overcome for the smaller signal levels expected from mini- 
halos. Note that a period of early reionization would wipe 
out the mini-halo population by raising the baryonic Jeans 
mass (Cen 2002). Also, early star formation would provide 
trace metal contamination, and ii Z ^ 10~^Zq (Bromm et al 
2001), the gas in mini-halos would be able to cool via metal 
line cooling and form stars. If the mini-halo signal is de- 
tectable in spectral pencil-beam measurements, it is prefer- 
able to have as small a bandwidth as possible (~ O.lMHz) 
to maximize intensity fiuctations due to the varying number 
of halos within a bandwidth along the line of sight. 
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5 CONCLUSIONS 

We have considered a number of foregrounds for 21cm emis- 
sion studies. We perform a robust minimal estimate of tem- 
perature fluctuations across the sky due to free-free emission 
by ionizing sources assuming only; (i) a number-weighted 
clustering bias with a cutoff for halos with T„ir < lO^K, 
(ii) a star formation history normalized to the metallicity 
of the z=3 Lyof forest, and (iii) an escape fraction of ioniz- 
ing photons into the IGM of no more than ~ few x 0.1. We 
conclude that the clustering component of temperature fluc- 
tuations will swamp the expected angular brightness tem- 
perature fluctations due to 21cm emission during the cos- 
mological Dark Ages by at least an order of magnitude. 
Performing 21cm tomography by detecting spectral features 
along a single line of sight requires developing scaled beams 
which sample exactly the same patch of sky at all frequen- 
cies. Otherwise, excess variance in foreground sampling due 
to the frequency dependence of the beam will swamp the 
signal. We also discuss radio recombination lines as possible 
contaminants. It is unlikely that this will be a problem for 
surveys outside of the Galactic plane, although this is still 
fairly uncertain. 

Detection of a 21cm spectral signature in absorption 
against a high-redshift radio loud source (Carilli, Gnedin & 
Owen 2002; Furlanetto & Loeb 2002) will not suffer from 
the foregrounds discussed here, except for the rather spec- 
ulative possibility that intervening high-redshift disks could 
give rise to a plethora of obscuring radio recombination lines. 
Of course, this 'foreground' would itself be a marvellous win- 
dow on the high-redshift universe! The main difficulty of 
absorption studies, of course, is whether such bright radio 
sources exist at early times, when structure formation is in 
its infancy. 
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